^1 VJI Wil c 



Sound propagation and force chains in granular materials 



Eli T. Owens and Karen E. Daniels^''^ 

Department of Physics, NC State University, Raleigh, NC, USA 27695 



o 

<N 

o 

Q 



o 

a 
s 

O 

> 

00 
O 

o\ 
m 

o 
o 



'i 



PACS 4 3.2 5. -X- 
PACS 4 5.7 0. -n- 
PACS 4 6.4 0.Cd- 



- nonlinear acoustics 

- granular systems 

- mechanical wave propagation 



Abstract. - Granular materials are inherently heterogeneous, leading to challenges in formulating accurate 
models of sound propagation. In order to quantify acoustic responses in space and time, we perform exper- 
iments in a photoelastic granular material in which the internal stress pattern (in the form of force chains) 
is visible. We utilize two complementary methods, high-speed imaging and piezoelectric transduction, to 
provide particle- scale measurements of both the amplitude and speed of an acoustic wave in the near-field 
regime. We observe that the wave amplitude is on average largest within particles experiencing the largest 
forces, particularly in those chains radiating away from the source, with the force-dependence of this am- 
plitude in qualitative agreement with a simple Hertzian-like model of particle contact area. In addition, we 
are able to directly observe rare transient force chains formed by the opening and closing of contacts during 
propagation. The speed of the leading edge of the pulse is in quantitative agreement with predictions for one- 
dimensional chains, while the slower speed of the peak response suggests that it contains waves which have 
travelled over multiple paths even within just this near-field region. These effects highlight the importance of 
particle- scale behaviors in determining the acoustical properties of granular materials. 



Introduction. - Sound propagation in granular materials 
differs from propagation in ordinary elastic materials in that 
there is a poor separation of length scales, particularly mani- 
fest in the branching networks of force chains which transmit 
stresses between particles. Continuum models \^^^, includ- 
ing effective medium theory (EMT), have failed to quantita- 
tively describe important features such as the dependence of 
the sound speed on pressure ||4j|6|. Particle-scale changes in 
the coordination number and force chains are likely responsi- 
ble for important deviations from these models, and have been 
the focus of numerical simulations on model amorphous sys- 
tems near jamming, where soft modes are important |7-9|. 

Hertzian contact theory underlines models of sound propa- 
gation, whether as the interaction potential between particles in 
discrete element simulations |T0| or in the calculation of an ef- 
fective bulk modulus in EMT |p]|. The contact force / between 
two particles is given by an equation of the form f oc S^, where 
S is the distance each particle is compressed and the exponent (3 
depends on the particle geometry. Two common idealized sit- 
uations are /3 = 3/2 for spheres and /3 = 1 for cylinders. The 
contact area a between the particles, which in EMT are perma- 
nently bonded, has a force dependence given by a ex / ^ , for a 
circular contact area in the case of spheres, and a ex / 2^ for a 
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rectangular contact area in the case of cylinders. The sound ve- 
locity c is governed by the stiffness s = df/dS of the contacts; 

for Hertzian contacts, this is 5 ex ^^ ex / ^ -In experiments 
on one-dimensional chains of identical spheres, |11| observed 
quantitative agreement with the expected c ex ^/s for a variety 
of materials. 

However, due to the force chain network present within gran- 
ular systems, the contact force, area, and stiffness are in gen- 
eral heterogeneous quantities, and measured acoustic signals 
arise from a superposition of signals traveling within this het- 
erogeneous medium. Furthermore, contacts between particles 
may open and close during the propagation 1 12"|, in contrast 
with the well-bonded assumption of EMT. Nonetheless, mea- 
surements of sound speed c as a function of system pressure P 
have observed the expected c (x P^^^ for 2D and 3D systems of 
spheres above a critical P for which non-Hertzian effects such 
particle asperities and contact rearrangements are minimized 
||2J|4|[T3|. However, nonlinear and non-affine effects prevent 
EMT from providing a complete treatment of the dynamics 1 6 1 . 

An additional complication lies in the observation that c 
can take different values within the same system, depending 
on whether the speed measurement is taken for time-of-flight, 
peak response, or harmonic excitation 1 24 1 . From this property 
of c, as well as the observation that a minute expansion of a 
single particle can cause large shifts in the acoustic response. 
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(b) t = 0.6 ms 




Fig. 1: (a) Image of photoelastic aggregate viewed through a pair of circular polarizers, 
region shown in (a), with white representing the areas with largest change in intensity. 
|http://nile.physics.ncsu.edu/pub/movies/gransound7| 



(b-d) Representative |A/| measured from 
A movie of the dynamics is available at 



|[T4^,^2?| proposed that the underlying force chain network is at 
the heart of the failures of EMT. Simulations by 1 16] directly 
addressed the heterogeneity and local dynamics using discrete 
element simulations of soft spheres. In a 2D system of ideal- 
ized Hertzian disks (/3 = 1), they observed that the propagating 
wave was in fact coherent and thereby insensitive to the force 
chain network. They interpret this result as arising from the am- 
plitude and speed both being independent of the local contact 
forces, as might be the case if the stiffness of the contacts were 
constant throughout the system for the special case of /3 = 1. 
Reconciling the failures of EMT with this unexpected insensi- 
tivity to the force chain network requires measurements of the 
propagation dynamics in real systems at length scales compa- 
rable to the particles and force chains. 

In order to address the particular role played by the 
inter-particle contacts and forces, we make spatiotemporally- 
resolved measurements in an experimental system in which 
changes in the force chain network are rendered visible. Our 
apparatus is filled with a two-dimensional (2D) aggregate of 
photoelastic |[T7]|T8| disks, allowing direct visualization of the 
sound path via polarizing filters and a high-speed camera. In 
addition, piezoelectric sensors embedded in a subset of the 
particles provide the higher temporal resolution necessary for 
speed measurements. This combination of techniques allows 
us to directly measure force, sound amplitude, and sound speed 
on the particle scale. We observe that the force chain network 
plays a key role in sound propagation, with the largest ampli- 
tude sound propagating along the strongest force chains. We 
find that a simple model of increasing contact area with inter- 
particle force provides semi-quantitative agreement with the 
observed trend. In addition, we find the sound speed of the 
leading edge of the signal in our 2D system is consistent with 
the speed found along a ID chain. However, the 2D speed of 
the peak response is not observed to be force dependent, and is 
slower than the corresponding ID speed. We use this result to 
explain the speed difference observed for 2D and 3D systems 
(241 by postulating that the leading edge of the signal travels 
only along the stiffest (fastest) path, while the peak response is 
a combination of sound that traveled along multiple paths and 
would therefore be independent of the force and slower than 
the leading edge of the signal. 



Experiment. - Our apparatus consists of a (29 x 45) cm^ 
vertical aggregate of (9(10^) photoelastic disks confined in a 
single layer between two sheets of Plexiglass (see fig. [T^). 
A pair of left and right circularly polarizing filters on oppo- 
sites sides form a polariscope, so that particles with larger 
force appear brighter and with more fringes. The particles are 
cut from Vishay PhotoStress material PSM-4, with diameters 
di = 9 mm and d2 = II mm in equal concentrations (d = 
10 mm); the zero-frequency Young's modulus of the material 
is £^0 = 4 MPa and the density is p = 1.06 g/cm^. Because the 
particles are viscoelastic, they have a stress-dependent dynamic 
modulus E{uj)\ for the 750 Hz signals used here, £^ ^ 50 to 
100 MPa 1 19|. The upper surface of the aggregate is free and 
we perform measurements near the bottom, where the Janssen 
effect |20||2T| eliminates vertical gradients in the average pres- 
sure. The average pressure in this region is 10~^£^o, with 
strong force chains typically 10~^£^o, and the system has an 
average packing fraction of 0.84 ± 0.01 To improve statistics, 
we perform experiments in 19 different particle configurations. 

A voice coil affixed to the bottom wall of the system sinu- 
soidally drives a flat platform of width 2J, providing a point- 
like source of acoustic waves. The driver is mechanically iso- 
lated from the apparatus to eliminate transmission of signals 
through the boundaries. The source pulse consists of 5 con- 
secutive sine waves with a frequency of 750 Hz, which corre- 
sponds to a wavelength A = 10 J to 20 J (based on sound speed 
measurement presented below). Prior to each experiment, we 
send an annealing sequence of 30 pulses so that the system set- 
tles into a state for which we observe repeatable measurements. 

We observe the wave propagation through the aggregate us- 
ing two techniques: the photoelastic response via a digital 
high speed camera, and the electrical response of particle- scale 
piezoelectric sensors embedded in a subset of particles (see 
fig [2J3). The former provides data at all particle locations, but 
only 4kHz temporal resolution; the latter provides data at only 
12 particles but 100 kHz temporal resolution. In addition, a 
3k X 2k digital camera provides higher spatial resolution im- 
ages for measurement of particle positions and forces. 

We image the photoelastic response of the aggregate using 
a Phantom V5.2 camera operating at 4 kHz with 512 x 512 
resolution. Since the change in the photoelastic signal dur- 
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Fig. 2: Sample signals showing correspondence between two mea- 
surement techniques, (a) Driver current (bottom), A/(t) for a repre- 
sentative photoelastic particle (middle) and piezoelectric voltage V{t) 
(top), with all values normalized to same scale. Gray bars illus- 
trate the intervals Atr and Atp. (b) Decay in average signal ampli- 
tude with distance from driver: AAi{r) for the photoelastic signal 
(O) and Av(r) for the piezoelectric signal (■). Solid lines show 
A oc exp~^/^ with ^ = 2.3 cm. Inset: particle with embedded piezo- 
electric sensor 



ing the pulse is weak, maximally ±10 out of 256 levels, we 
transmit the 5 -wave pulse 50 times at 20 ms intervals to al- 
low the previous pulse to dissipate. The frame-by-frame aver- 
age intensity I{x,y;t) of these 50 pulses is used to measure 
spatially-resolved dynamics. In each frame, we determine the 
locations where compression is present by examining differ- 
ences in the image brightness with respect to the initial frame, 
AI{x,y]t) = I{x^y]ti) — I{x^y;to), as shown in fig. [TJ)- 
d. The field Aai{x^ y; t) records the maximum amplitude ob- 
served at each particle location. Note that for a measured speed 
of c = 200 m/s, the 4 kHz frame rate averages the photoelastic 
response over a distance of approximately bd in the direction 
of propagation. 

In addition, we utilize embedded piezoelectric sensors to 
provide higher temporal resolution for a subset of the parti- 
cles. Piezoelectrics were chosen because they can fit within 
a single particle, produce a voltage proportional to the applied 
stress, and have a frequency response well above the frequen- 
cies we measure. Our sensors are 2 mm thick lead zirconate 
titanate piezoelectric ceramics from Piezo Systems, Inc., cut 
to 4 mm X 4 mm squares; this adds 30% to each particle's 
mass. The configuration of the 12 sensor-particles is shown in 
fig.[T^, with the first row 5 cm from the driver and the remain- 
ing rows spaced 3.5 cm apart. Since our piezoelectric sensors 
are only sensitive to stresses applied across their thickness, it is 
important to orient them so that their face is perpendicular to 
a line connecting the center of the driver with the center of the 
sensor. The voltage measured on any particles not in this ori- 
entation is corrected with the appropriate cosine factor during 
post-processing. 

During each train of pulses, we record the time-varying volt- 
age produced by each of the twelve piezoelectric sensors (see 
fig. [2^). Due to the high impedance of the piezoelectric sensors, 
an op-amp buffer and a 10 M£t resistor are connected in parallel 
with the piezoelectric to lower the impedance of the circuit and 
reduce signal crosstalk. Three quantities are measured for each 
sensor and the driver: amplitude Ay, rise time tr, and the mid- 



dle peak time tp. We measure amplitude Ay by fitting a sine 
function to the middle three peaks in V{t). We measure tr as 
the time at which each signal rises 15% above its resting value, 
and tp as the time at which the maximum value of the middle 
peak is attained for each signal. The intervals Atr and Atp are 
measured for each sensor with respect to the driver signal, and 
are used to calculate the speed of sound c along the observed 
sound path. 

The two amplitude measurements, photoelastic Aa/(^, y] t) 
and piezoelectric Ay{x^y]t) provide complementary informa- 
tion: the former with better resolution in (x, y) and the latter 
with better resolution in t. As can be seen in fig. [2^, A/(t) 
for a single particle shows the expected 5 -peaked signal sent 
by the driver and measured at a nearby sensored particle. In 
fig.[2j), for all 19 configurations, we average Ay(r) over the 3 
sensors in each row and azimuthally average A/^i{r). We find 
that {Ay{r)) and (AA/(r)) exhibit the same exponential de- 
cay. As such, the two amplitudes provide the same information. 
Note that this exponential decay is due to the dissipative nature 
of both the viscoelastic photoelastic polymer and the frictional 
granular contacts, and suppresses the multiply-scattered (coda) 
response ||14[|22|. 

These particle-scale measurements of sound amplitude and 
propagation speed can be correlated with the measured particle- 
scale forces in order to quantify the effect of the force network 
on the propagation dynamics. In order to find the force F on 
each particle, we locate all particles on a high-resolution im- 
age and calculate the squared average intensity gradient | V/p 
within each particle. This empirical method 1 18 1 is calibrated to 
known forces in a separate test cell. Three calibration regimes 
were used for low-, mid-, and high-force particles. 

Results. - 

Contact Force Law. We directly measure the contact force 
law for our particles by compressing a single particle between 
two plates with the force applied along the particle's diameter, 
fig. |4^ shows the measured force for a range of applied dis- 
placements. Within a range which corresponds to the typical 
forces in the full experimental system, we find that the contact 
force law is approximately of the form 



/ oc ^^/^ 



(1) 



While this deviates from the Hertzian prediction that co-planar 
cylinders follow a linear contact law, we note that slight out- 
of-plane rotations of the particles lead to imperfect alignment 
at their contacts. The value /3 = 5/4 can therefore be inter- 
preted as lying between (3 = 1 (for rectangular contacts be- 
tween co-planar cylinders) and /3 = 3/2 (for circular contacts 
between perpendicular cylinders) (23). As will be seen below, 
this deviation from a linear contact law has several detectable 
consequences. 

Linearity and nonlinearity. The photoelastic particles al- 
low us to directly visualize an important particle-scale nonlin- 
earity present in the system. Due to the non-zero Poisson ratio 
of the polymer material {y ^ 0.5), a particle under dynamic 
compression expands in the perpendicular direction, potentially 
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Fig. 3: (a) Detail of force chains approximately 8 J from the driver, 
with particles 1 and 2 exhibiting weak contacts, (b) | A/| during sound 
propagation for the same region. New contacts to particle 1 (►) and 
the new chain through particle 2 open and close repeatedly during each 
wave cycle. 
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Fig. 4: (a) Measured contact force law compared with / oc ^^^^ and 
f (X 6 (solid lines), (b) {Ay) found from all twelve piezoelectric 
sensors in a single pack as a function Ain. Ay and Ain are in arbitrary 
units of pressure. Error bars are standard error and the solid line shows 

Ay OC Ain. 



creating new (but transient) contacts during the transmission of 
a wave. The photoelastic particles render the creation of such 
new mechanical contacts visible, as can be seen in fig. [3] In 
the static image of the force chains (a), particles 1 and 2 are 
not part of the strong force chain network. However, during 
the transmission of the pulse (b) each of these marked parti- 
cles exhibits a measurable |A/| before returning to its initial, 
un-stressed state. This observation that particle contacts can 
change during sound propagation violates the well-bonded as- 
sumption of effective medium theory 1 1 1 because the connec- 
tivity of the network transiently changes due to the passage of 
the wave. Similar, but permanent, changes to the force network 
were presumably at work in experiments where minute thermal 
expansion |24| or consolidation of loose packs |25 1 modified 
the signal coda. 

In spite of the inherent nonlinearity of such events, we 
nonetheless observe an approximately linear response in the 
average propagation amplitude of the system. As shown in 
fig. [4]3, we measure the input amplitude Ain as the pressure 
generated by the driver as it compresses the particles, and the 
output amplitude as the average signal {Ay) received by all 12 
piezoelectric sensors in a particular configuration. Since the 
piston is driven sinusoidally, it has a maximum displacement 
of J = ^^^, where Fnou is the electrical force from the coil 
(measured via the current), uo is the driving frequency, and m is 
the mass of the piston plus some effective mass of the system. 
For a given uo, system configuration, and Fcoiu this unknown 
mass remains constant. Therefore, 5 oc FcoU is not only the 
displacement of the piston, but also the compression applied to 
the first particle. Using this relation to provide 6 in eq. [T] we 

5/4 ^"^ 

find that the applied pressure Ain ^ ^coU • ^^ ^^ increase the 
current to the coil, we can thereby measure {Ay) as a func- 
tion of Ain and find an approximately linear relationship. This 
suggests that at the driving magnitudes investigated here, the 
transient contacts are not a large enough effect to detectably 
destroy an approximately linear response. 

Amplitude. As can be seen in fig. [T]^-d, the dynamic re- 
sponse is quite similar in structure to the original force chain 
structure shown in fig. [T^. To quantify this effect, we measure 



both the local force F on each particle and the maximum am- 
plitude Aa/ of the wave passing through that particle. For all 
particles in each of the 19 configurations, we bin the particles 
by F and calculate the average (Aa/), where A is corrected 
for the exponential decay shown in fig. [2J3. As can be seen in 
fig. [5] the particles with stronger forces on average experience a 
larger amplitude, which corresponds to the visual observation. 
We can understand this amplitude trend in terms of the mea- 
sured contact law in eq. [T] The Hertzian contact area a formed 
between two such disks pushed together by a pair of oppos- 
ing forces / is a oc P^^. If the contacts between particles are 
the main conduits for sound propagation, then we would also 
expect the magnitude of the sound to vary as Aa/ cx F^/^. 
This relationship is consistent with the data shown in fig. [5] 
particularly outside the low-force regime. This contact model 
is an over- simplification in several important ways. While the 
Hertzian contact theory was developed for diametric loading, 
the multiple contacts on each particles arise at various angles. 
Furthermore, while the sound propagates radially outward from 
the driver, the force chains form arcs which are not in general 
aligned with this direction (see fig.[l]). As such, a more sophis- 
ticated treatment would consider only the radial component of 
the stress tensor. The fact that we recover an approximate scal- 
ing of the amplitude points to the robustness of the trend. 

Sound Speed. We measure sound speed via the signal 
travel time from the driver to the piezoelectric sensors (see 
fig. [2^) along the dominant propagation path from the driver 
to each sensor from the high-speed movies (see fig.[T]3-d). As 
in | p^[24| , we record both the travel time for the initial rise 
(At^) and for the peak response (Atp). The distance travelled 
is taken to be the observed path (rather than Cartesian distance) 
since this is the path over which we observed the largest ampli- 
tude sound propagates. This path is selected by hand from the 
movies; we measure the path length L and calculate the two 
sound speeds, c^ = L/Atr and Cp = L/Atp. In addition, we 
measure the average force F from the average force within the 
particles along the observed path. Due to the exponential decay 
in the amplitude of the signals, quantitative speed data is only 
available for the inner two rows of sensors located at r < A. 

In fig. l6k, we observe that Cr is on average faster than Cp 
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Fig. 5: Average scaled sound amplitude (Aai) as a function of par- 
ticle force F, with each data point averaged over several hundred 
measurements. Symbols (Q, V, D) are for calibrations in the low-, 
mid-, and high-force regimes, respectively, with error bars for indi- 
vidual measurements in each regime shown at bottom of plot. Line is 
A oc F^/^ 
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Fig. 6: (a) Sound speed measured by the piezoelectric sensors as a 
function of the average force F along the observed path. Large sym- 
bols are for sensors at r = 5 cm, small symbols at r = 8.5 cm. The 
speed of the initial rise Cr (gray •) and of the peak Cp (gray A) with lin- 
ear fits to each. Black symbols are averages of the same data, binned 
by force. Linear fits are: Cr{F) = (106 it 46) m/s +(130 ± 97) s/kg 
F Cp{F) = (93 ± 38) m/s -(20 ± 80) s/kg F with a confidence in- 
terval of 95%. (b) Comparison of the ID (gray) and 2D (black, copied 
from (a)) speeds, Cr (•) and Cp (A), with c oc /^^^° for compari- 
son. Bulk speeds ( x ) are calculated from measured dynamic modulus 
under fixed compression 1 19|. 



(as was previously observed in 1 16 24j, see Discussion for the 
significance of this result), Cr and Cp are both independent of 
distance from the driver as expected, and that there is signifi- 
cant scatter in the data which impedes our interpretation of the 
force-dependence. As most of the scatter in fig. [6^ is due to 
uncertainty in the photoelastic force measurement, we bin the 
data by force in order to clarify the force-dependence, and com- 
pare it with both Hertzian predictions and speed measurements 
in one-dimensional (ID) chains. 

We performed sound speed measurements on ID chains of 
the 1 1 mm particles confined by a known force within a linear 
channel centered above the driver. For Hertzian contacts, we 
expect c (X ^/s ex /^^ ex /^/^^, which provides reasonable 
agreement with the measured sound speeds, Cr and Cp, on the 
ID chain shown in fig.jSJ). Notably, these speeds are faster than 
the bulk speed predicted by the zero-frequency Young's modu- 
lus (c = ^/Eo/p = 61 m/s). However, E depends on both the 
driving frequency and degree of pre-loaded stress due to vis- 
coelastic stiffening. Independent measurements of the dynamic 
E at typical values of uj and / (T9J, predict a speed consistent 
with what we observe. 

The agreement of the ID and 2D sound measurements of 
Cr (see fig. [6|)) suggests that the Hertzian contact model pro- 
vides a good description of the weak force-dependence. This 
observation contrasts with numerical simulations using ideal 
(13 = 1) Hertzian disks and Coulombic friction, where |T6| ob- 
served no such force-dependence in either speed or amplitude 
measurements. We attribute this difference to our measured 
force law (eq. [T]) which has /3 = 5/4; in a 3D aggregate of 
Hertzian spheres (/3 = 3/2), similar force-dependence might 
be expected. Additionally, we note that in 2D Cp is slower than 
Cp in ID and is approximately independent of the force (see 
Discussion below.) 

Discussion. - These experiments are performed on a 
highly dissipative viscoelastic material within a near-field 
regime (r <C A). Nonetheless, we are able to reproduce features 
seen in more conventional 3D granular materials (c^ > Cp) and 
ID chains (c ex ^/s) which indicate a degree of universality 
among systems of quite different dimensionality, particle ge- 
ometry, and material properties. In part, this is due to the fact 
that the viscoelasticity only affects the constant of proportion- 
ality in / ex ^^, while /3 depends only on the contact geometry. 
However, it is possible that such additional stiffening enhances 
the importance of the force chains. In spite of these limitations, 
the use of these photoelastic materials permits us to address the 
importance of particle-scale heterogeneities to sound propaga- 
tion. The application of these techniques outside the near-field 
regime will require the use of less-dissipative materials. 

As shown in fig. |5j we observe that the local sound ampli- 
tude increases with the force on the particle, with a shape that 
suggests the the effect is due to increasing contact area (A ex 
a ex /^/^). This finding is in contrast with a lack of amplitude- 
dependence in simulations of idealized (f3 = 1) disks. 1 16|. For 
realistic deformable particles (rather than merely overlap with 
a repulsive potential), the contact area should still increase with 
force: a ex / 2^ ex /^/^. 
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In experiments on sand p4| and in numerical simulations 
p6| , the measured value of c depended on the method: the 
speed at the leading edge of the signal, c^, was observed to be 
several times faster than that of the peak response, Cp. By def- 
inition, Cr represents the arrival of the first signal, that which 
arrived along the fastest (stiffest) path. Through the use of pho- 
toelastic particles to measure the dominant sound path, we ob- 
serve that these fast/stiff paths correspond to force chains. Be- 
cause Cr is sensitive to the particular force chain configuration 
through s (X f ^ , it shows evidence of (weak) force depen- 
dence (see fig. [6})). In contrast, the peak response measured by 
Cp is composed of signals which have travelled via many dif- 
ferent paths of varying stiffness and amplitude. Because this 
portion of the wave has traveled on many different paths, the 
speed is on average independent of the force along the fastest 
path, as seen in fig. [6^. This may be related to simulations 
which show Cp saturating to a maximum value as a function 
of r f26l. Also note that Cp is faster in ID chains than in the 
2D system, since all signals must travel on the same particle- 
scale path. The remaining speed deficit is thereby due to either 
dispersion or path-effects at the sub-particle scale. 

Thus, prior observations of c^ > Cp are likely due to c^ mea- 
suring sound that travelled along the (stiffer) force chain net- 
work. A consequence of the force-chain sensitivity of c^(/) 
is that the leading edge of the wave front will not form a co- 
herent arc propagating outward, but will instead be irregular. 
The behavior along each segment of force chain may be rem- 
iniscent of sound propagation in ID granular chains |11|, but 
with a network of connections between the chains. Interest- 
ingly, in the limit of zero confining pressure (no force chains), 
the observation that the intercepts of c^(/) and Cp{f) are ap- 
proximately equal would appear to suggest that the compres- 
sional wave speed c takes a single value in this limit. However, 
in studies of granular acoustics near this jamming transition, 
additional nonlinear effects are present which are relevant to 
understanding the scaling of shear and compressional waves 

I2L28J- 

Conclusions. - In conclusion, we observe that both the 
amplitude and speed of sound are influenced by the force 
chains, and that these chains may change due to the transmis- 
sion of the sound pulse. This highlights the role of force chains 
as more than just a contact network between grains. Local 
properties of the material are crucial to local sound propaga- 
tion dynamics, and mean field models therefore are unlikely 
to capture the full behavior of the system. The dependence of 
these results on the wavelength of the disturbance | [29| , either 
above or below the length scale set by the force chains, may 
provide a means to acoustically probe the force chain network 
in non-photoelastic 3D systems. 
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